TAN-FNT-96/003 



hep-th/96l2m 



Covariant Quantization of the Skyrmion 



M. Kruczenski 

Departamento de Fisica, Comision Nacional de Energia Atomica, Av. del 
Libertador 8250, 14-29 Buenos Aires, Argentina. 



and 



L.E. Oxmani 



Instituto de Fisica, Universidade Federal do Rio de Janeiro, CP. 68528, Rio 
de Janeiro, Brasil. 



Abstract: We obtain the internal degrees of freedom of the skyrmion (spin 
and isospin) within a manifestly Lorentz covariant quantization framework 
based on defining Green functions for skyrmions and then, the S-matrix via 
LSZ reduction. Our method follows Frohlich and Marchetti's definition of Eu- 
clidean soliton Green functions, supplemented with a careful treatment of the 
boundary conditions around the singularities. The covariant two-point function 
obtained propagates a tower of spin equal to isospin particles. Our treatment 
contains the usual method of collective coordinates, as a non-relativistic limit 
and, because of the new topology introduced, it leads, in a natural way, to the 
inequivalent (boson/fermion) quantizations of the SU(2) skyrmion. 



1 Introduction 



Since the work of Witten |T]J^, the Skyrme model ^ has received renewed 
attention as the low energy limit of QCD, as it is supposed to describe QCD 

^ E-mail: martink@tandar.cnea.edu.ar 
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in the large number of colors (Nc) limit. Recently, B. Moussallam gave 
further support to the soliton picture showing that the low energy expansion 
of chiral perturbation theory supplemented with a large A^^^ analysis is 
compatible with a stable soliton. 

The large Nc limit can be thought of as a semiclassical expansion, as Nc ap- 
pears in the combination h/Nc 0. One starts from a minimum of the potential 
which is the classical soliton profile. The soliton solution has six zero modes 
corresponding to translations and rotations (or isorotations). Quantization of 
these zero modes gives rise to linear momentum, spin and isospin. An analysis 
of the wave functions reveals the existence of a tower of particles with spin 
equal to isospin. 

Also, as the mass of the skyrmion scales as Nc, in the large Nc limit a non- 
relativistic treatment of the skyrmion is justified. Although this is correct in 
principle, it is not always desirable to break Lorentz invariance as one looses 
the possibility of applying standard knowledge from quantum field theory. 
For example, in reference 0, the equivalence to a relativistic field theory 
was advocated as an important point to solve the so called Yukawa problem. 
Therefore it is useful to introduce a manifestly Lorentz covariant quantization, 
defining Green functions for skyrmion fields. We believe that it is conceptually 
better to have a definition for these Green functions and, a posteriori, to take 
the non- relativistic limit. 

Frohlich and Marchetti p| introduced Euclidean Green functions for skyrmi- 
ons following previous developments on soliton quantization P,p!0|. In refer- 
ence 



11] similar ideas were considered in the context of gauge theories. 



In the present paper we will consider the covariant approach of ref. [§] but, in 
order to obtain the skyrmion internal degrees of freedom (spin and isospin), 
we will prescribe the boundary conditions which the field satisfies around the 
singularities. This treatment will allow us to relate the covariant quantization 
approach to more traditional ones for skyrmion quantization^. Also, this co- 
variant quantization is a natural candidate for the skyrmion quantum field 
theory underlying the work of N. Dorey, J. Hughes and M.P. Mattis 

In section |^ the Green functions defined in are reviewed. In section ^ the 
boundary conditions are discussed, showing that they bring in new topology. 
In section ^ the propagator is shown to transform correctly under Lorentz and 
isospin transformations. Further arguments are given in section |^ to show that 
the propagator describes a tower of spin equal to isospin particles. In section | 
we give the LSZ rules to construct the S-matrix and discuss the relation of 
our calculation to the computation of ref. . In section ^ our conclusions are 
given. Finally, in an appendix, we summarize the definition of arbitrary spin 

^ For a review see 
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propagators due to Weinberg 
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2 Green functions 



The Skyrme model is defined by the following Lagrangian 



f 1 
L = ^TidM^df'U + -Ti { U^dM, U^dM 



where [/ is a matrix belonging to SU(2). The lagrangian is invariant under 
global SU(2)lxSU(2)j^ chiral symmetry transformations U QlUqIi, Ql, Qr G 
SU(2). This symmetry is spontaneously broken. If one chooses the vacuum 
expectation value of U to be (U) = 1, the surviving symmetry is the diagonal 
subgroup SU(2)y {gi = gn). A convenient parametrization of U is 

[/ = exp(i0ar"), (2) 

where are the Pauli matrices, and the fields 0a describe the pions which have 
zero vacuum expectation value {{U) = 1 ^ {(pa) = 0). One can also consider 
a pion mass term and higher derivative terms such as —e^B^^B^, where ee is 
a constant and B^ is defined in eq. (|^) below. From now on, we will work in 
Euclidean space and will use upper indices (/x = 1, 4) to denote space-time 
components. In the following, 4-dimensional Euclidean rotations will be called 
Lorentz transformations for the sake of brevity and to distinguish them from 
3-dimensional rotations. However this is an abuse of language and the reader 
should keep in mind that we are working in Euclidean space-time. 

Static configurations with finite energy must satisfy f/ — 1 at infinity, this 
amounts to compactify space to S'^. Thus, a static configuration defines a 
map — > SU(2), that is, an element of 113(5'^) = Z. The baryonic number 
B, defined as the integral of the topological current 

B = j B\ (3) 



is the winding number associated with 113(5''^) and therefore is conserved. 
Configurations with B = 1 describe nucleon states. The energy of the B = 1 
static configuration that minimizes the action is concentrated in space; then, 
this configuration can be approximately treated as a rigid body rotating in 
space and also in internal space. Spin and isospin arise quantifying these col- 
lective degrees of freedom M . 
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If there are quantum states associated with the classical soliton then it should 
be possible to define Green functions such as 

^n,...n„ai...a,.(^^ • • • X^, ^/i . . . y,) = (T (rH^l) • • • {Xp)<PaAyi) ■ ■ ■ KiVr))) . 

(4) 

where ip"'^{xi) (formally) represents an operator destroying a soliton of winding 
number nj at Xi or creating an antisoliton — at Xi. 

Following P] we introduce a path integral over fields defined on space-time 
compactified to S"^ and with the points xi...Xp removed (Mi p = {R'^ U 
{oo})\{a;i, . . . ,Xp}). The manifold Mi p is contractible to a bundle of p — 1 
spheres with one common point, implying that maps from Mi p to SU(2) 
are classified by p — 1 independent integers. These integers can be written in 
terms of p integers Ui, satisfying rii + ■ ■ ■ + Up = 0, defined by 

ni = JdS n'^B'^ ,i = 1...P, (5) 

where Si is a three sphere enclosing Xi which does not enclose any other point 
Xj (j ^ i), is the normal to Si and dS is the surface element on Si. Due to 
the fact that U is constant at infinity, we have rii + ■ ■ ■ + = 0. 

The skyrmion Green functions in eq. (^ were defined by Frohlich and Marche- 
tti in ref. |^. Including the pion field definition (H), the mixed pion-skyrmion 
Green functions read 

^n,...n„a,..a,.(^^^ • • • , Xp, . . . y,) = j P[f/]</)a, (Z/l) • • • 4>aM ^M' J '^P]) , 

M 

(6) 

(see 1^ for an interpretation in terms of open line defects). U = exp{i(j)^T^) 
takes the value [/ = 1 at space-time infinity and has topological number 
around Xi {i = l,...,p). A rigorous definition of the path integral would 
involve first defining it on a lattice and then letting the lattice approach the 
continuum . Alternatively we can define the path integral by its perturbative 
expansion around a saddle point. In both finite ultraviolet cut-off A 

will be required since the Skyrme model is non-renormalizable. In the case of 
chiral perturbation theory the limit A —>■ oo can be taken by absorbing the 
infinities in the higher order derivative terms, which are neglected at large 
distances. 

In order to clarify the meaning of eq. (^ let us consider the soliton propagator 

G-i'+i(xi,X2) = J V[U]exp{- J C[U]). (7) 

A/l2 

The path integral is performed over fields U which have winding number +1 
around X2 and —1 around Xi. The situation is depicted in figure 
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Fig. 1. To define the propagator, configurations over R'^\{xi,X2\ are considered. 
See the discussion in the text. 

Consider a spatial section at constant time t = ta- By closing this section 
at t = —oo we see that it has topological number 0, that is, there are no 
topological excitations at t = ta- The same argument for a section at t = tb 
shows that it has winding number +1 and, at t = t^, the winding number is 
again 0. Then, we see that a. B = 1 soliton is created at X2 and destroyed 
at Xi. Furthermore, if we look at the time inverted process, we find that 
section t = tf, has topological number —1, because definition (|) picks up 
a minus sign when the normal is reversed. Thus, the figure also describes an 
antisoliton propagating backwards in time, that is, at X2 a soliton is created or 
equivalently an antisoliton destroyed, the effect expected from a quantum field 
operator associated with the antisoliton. Remember that a field destroys the 
particle and creates the antiparticle whereas ip'^ does the opposite. 

The following step is to find a saddle point configuration for our path integral. 
This can be done proposing a (time-dependent) Skyrme ansatz 

t/ = exp(if(u,v)xar^), (8) 



where Xa = Xa/r, r = yxf + + X3, and u = ef-„r, v = ef-^t are adimen- 
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sional. Minimizing the action, gives an equation for f{u,v): 




^ (if ^ du"^ ^ 

+2w£ - sin(2/) ( 1 + ) = 0. (9) 

The boundary condition U = 1 at space-time infinity requires that f{u, f ) — > 
when V —>■ ±oo or m — »• cxo. At m = 0, for U to be single valued, / must be 
or TT. Continuity of / (on the manifold M12) and the requirement that U must 
have topological number 1 around X2 and —1 around xi leads to 

f{u,v) = if t> ^ ±00 or ti — >• cxD, 

/(0,f) = if t> < efT,X2 or V > efj^x^, (10) 

f{0,v) = TT if ef^xl < V < ef^xl- 

These boundary conditions respect the cylindrical symmetry of the ansatz. In 
the next section we will discuss a more general situation. 

The numerical solution to eqs. (||) and (|10D is depicted in figure | where we 
see that / (and so the fields (pa) is approximately zero before or after xf 
and is almost equal to the static soliton profile in between. This supports the 
picture that a skyrmion is created at X2 and propagates towards Xi where 
it is destroyed. Furthermore, since the solution is almost time independent 
in the time interval (a;|,a;f) and equal to the static profile, the action is ap- 
proximately given by 5* = Mci\x2 — Xi\, where M^i is the energy of the static 
solution. Hence, the saddle point approximation to the propagator contains a 
factor exp(— Mci|a;2 — a;i|) which is the exponential behavior of the Euclidean 
propagator for a particle with mass M^. In section 5 we will improve this 
approximation, obtaining the correct, spin equal to isospin, particle spectrum. 

Let us now consider spheres Sr of radius R < \x2 — xi\ around X2- All these 
spheres have topological number 1 but the topological density, which in this 
case is given by p(a;^) = -B'^(x — X2)^/\x — X2\, will be differently distributed. 
Replacing the ansatz of eq. (^) in the last expression it follows that the density 
is only a function of 6, the angle between x — X2 and the 4-axis, and is given 
by 

ZTT^ sm (t') 

The numerical solution f{u, v) as a function of 6 can be obtained, by interpo- 
lation, from the formula 

f{e) = f{Rsm{e),e^xt + Rcos{e)). (12) 

The result for several values of R (in units of (e/^)~^) is depicted in figure |^. 
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This figure is available at 

http:/ /www.tandar.cnea.edu.ar/preprints/FNT/96/003.tex 



Fig. 2. Numerical solution for f{u, v) that minimizes the action (u = ef-,^r, v = ef.„t). 
In the figure, \x2 — xi\ = 6(e/vr)~^. 




0.0 0.5 1.0 1.5 2.0 2.5 3.0 



e 

Fig. 3. / as a function of for different radius, corresponding (from the upper to 
the lower curve) to R = 0.5,1.0,1.5,2.0,2.5,3.0 in units of {ef-^)~^, in the case 
\X2 - Xl\ = 6(e/7r)"^ 

In this figure, we see that near the singularity {R 0) f{6) tends to vr — 
which gives a uniform topological density (cf. eq. (|n|)). As one gets further 
from the singularity the topological density concentrates around the North 
pole of the sphere giving rise to a soliton moving towards xi. This is again 
the graphical picture corresponding to a soliton appearing around X2. This 
behavior of the solution (unlocalization for small R, and localization for large 
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R) is in agreement with studies of the static ansatz on a sphere and is in 
correspondence with the fact that at short distances the chiral symmetry is 
unbroken whereas it is broken at large distances. 



We end this section noticing that the behavior f{6) -k — 9 near X2 is equiv- 
alent to 

U{x)^-Xl{x), X,{x) = y ^ + ii^^ra, (13) 

\X — X2\ |X — X2I 

where are the Pauli matrices. A similar analysis around Xi, reveals that U 
behaves as Xi, where Xi is defined as X2, but with X2 replaced by Xi. These 
behaviors will be useful in the next section, when discussing the boundary 
conditions to be imposed on the fields. 



3 The boundary conditions 



The numerical saddle point obtained in the previous section is not invariant 
under rotations or isorotations, that is, there are several other saddle points 
or equivalently, several classical trajectories. This is due to the fact that we 
have not yet specified the boundary conditions satisfied by the field around 
xi and X2. 

In order to do so, we extract two small balls Bf and i?| of radius e — > 
around Xi and X2 respectively. Then, we need to specify the values taken by 
the field U on the surfaces = dB\ and S'| = dB^ of these balls, satisfying 
the condition that U must have topological number +1 around X2 and — 1 
around Xi. 

In principle, any function U : SU(2) can be chosen, but we will restrict 

ourselves to a minimum set of boundary conditions which close under Lorentz 
and isospin transformations. At the end of section |^ we saw that near X2 the 
numerical solution behaves as — X], resulting in a uniform topological density. 
Performing Lorentz transformations we obtain, around the singularity, the 
behaviors Si{-X2ysl {S^ S2 G SU(2), see eq. (^) below) without altering 
the uniform topological density. This suggest that all the boundary conditions 
Si{—X2y Sl should be treated on an equal footing to obtain a representation 
of Lorentz and isospin sjTiimetries. Thus, the boundary conditions are taken 
to be 

[71^3 = AiXiBl 

U\s^ = A2{-X2yBl (14) 
where Ai, Bi, A2, B2 belong to SU(2) and Xi,X2 are the SU(2) matrices de- 
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fined by 

4 4 a a 
Xl,2{x) = -. '-r + ^Ta, (15) 

As an aside remark, note that, in terms of tlie vector field = U'^d^U , these 
conditions look like the insertion of a zero-radius instanton (anti-instanton) 
at xi (a;2).Q 

Each pair Aj, Bi determines a map SU(2). In fact, as exactly the same 

map is determined by —Ai, —Bi, the boundary condition is defined by a point 
in the manifold SU(2) x SU(2)/Z2. These boundary conditions represent the 
internal degrees of freedom of the skyrmion fields. 

Given two points pi and p2 on SU(2) x SU(2)/Z2 we define a propagator 

/ \ 

GixuPi;x2,P2) = \}ml V[U]exp - J C[U] + 2C, , (16) 

\ A/i2 J 

where U at Sf is given by pi and at 5*1 by p2, according to eq. (0) (at space- 
time infinity we take U = 1, the vacuum configuration). = — (Svr^/e^) log(e) 
is added to cancel the infinite contribution to the action due to the singularity 
and can be thought of as wave function renormalization of the skyrmion fields. 
Including higher order derivative terms in the action would produce worse sin- 
gularities and should be properly redefined; the inclusion of, say, the sixth 
order term —e^B^^B^ requires that = — (67r^/e^) log(e) — e^{2T('^t^)^^ . Of 
course, at short distances, the most significant terms are precisely those ne- 
glected at large distances, reflecting the fact that near the singularities the 
lagrangian in eq. (|Tp is not adequate. However this does not modify the sym- 
metry arguments and, far from the singularities, the classical solution is in- 
sensitive to the higher order terms. 

Now, we will show that the imposed boundary conditions lead to a natural 
discussion of the different possible quantizations of the Skyrme model. Up to 
now, as Mi2 = (-R^ U {oo})\{a;i, X2} is contractible to S^, the configurations 
U : M12 SU(2) are classified by an integer. However, when we fix the initial 
and final points we have a map 

f/ : 53 X [0,1] ^ SU(2), (17) 

because M12 can be deformed to a cylinder with base S^. Here, dB\ and 
9-B| are the bases of the cylinder and the parameter in [0, 1] represents the 
continuous deformation of (9i?| into dBf, passing through sections between xi 
and X2 (remember that space-time infinity is identified). Two such maps are 

^ In our case we do not have a gauge theory, but one can define the Green functions 
in terms of a field U coupled to a singular as described in HI 
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not necessarily homotopic. The situation is similar to that of a plane with a 
hole. Paths with free end points are all homotopic but paths with fixed end 
points are classified by ni(S'^) = Z . 



In our case suppose we give two configurations 



f/i,?72 : 5^ X [0,1] ^ SU(2) 



(18) 



which satisfy the same boundary conditions at and 1. They can be used to 
define a map V : S"^ ^ ^ SU(2) as: 



It is easy to verify that V is continuous at 1/2 and periodic ([/(X, 0) = 
[/(X, 1)), provided that U\ and satisfy the same boundary conditions at 
t = 0, 1. This map is the analog of the closed path which can be constructed 
from two paths with the same end points, in the simple example of the plane 
with a hole. 

Two configurations \J\ and t/2 are homotopic if the map f/ : 5'^ x 5^ ^ SU(2) 
can be extended to a map S*^ x ^ SU(2), where D2 is a two dimensional 
disc {pD^ = S^). In general, there is an obstruction to such an extension 
given by 114(811(2)) = Z2 0, this means that there are two different homotopy 
classes. A way of distinguishing them is to trivially embed SU(2) into SU(3) 
and evaluate the Wess-Zumino-Witten term. This gives for the trivial class 
and vr for the non-trivial one. If we were considering SU(3) solitons all the 
paths would be homotopic and the Wess-Zumino-Witten term would give a 
different contribution for each path 0. We would like to point out that the 
need of considering an extension of a map U : x ^ SU(2) to a map 
U : X D2 SU(2) to study the different possible inequivalent quantizations 
of the Skyrme model, which in was proposed by imposing periodic boundary 
conditions to calculate the partition function, here also arises naturally from 
the Frohlich and Marchetti's formalism supplemented with fixed boundary 
conditions around the singularities. 

The path integral defining the propagator can then be split into two inte- 
grals, one over each homotopy class. A phase equal to or vr can be assigned 
to the non-trivial class, which corresponds to adding or subtracting the two 
contributions: 



The plus (minus) sign will correspond to quantizing the soliton as a boson 
(fermion) (see section |^) . 




if t e [0,1/2] 
t)) if tG [1/2,1]. 



(19) 



G{xi,pi;x2,P2) = G^^\xi,pi;x2,P2) ± G^^H^^i, Pi; 3^2,^2) • 



(20) 
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Now, suppose that we have a path in some topological class which ends at 
a boundary condition associated to ^2,-82! if we adiabatically move this last 
boundary condition, leaving A^, Bi fixed, to get a boundary condition associ- 
ated to — A2, — -B2 (which represents the same point as ^42, B2 and —^2,-52 
are identified), the new path will be in a different topological class. This is 
due to the fact that the new path can be obtained by composing the original 
path with a "closed" path, based on (9-B|, that starts at (^2,-82) and ends at 
{—A21 —B2). This last path is non-trivial and so the original path will change 
from homotopy class. Thus, changing adiabatically the final condition from 
(^2,^2) to {-A2, -B2) interchanges G'(^) and G^^). 

This result allows us to continuously extend the boundary conditions from 
SU(2) X SU(2)/Z2 to SU(2) x SU(2), defining a new function G according to 

(5(^1, A, Si; X2, A2, B2) = G{xi, -A, -Si; X2, -B2) 

= G^^\xi,pi;x2,P2), 
G{xi, Ai, Bi; X2, -A2, -B2) = G{xi, -Ai, -Si; ^2, A2, B2) 

= G'(2)(a;i,pi;x2,P2), (21) 

where Ai^2-, Si,2 and —Ai^2, —Bi^2 both correspond to pi,2 in SU(2) x SU(2)/Z2 
and we take the definitions to be valid when Ai, ^2 > ^. In terms of G the 
original propagator can be written as 

G{xi,pi;x2,P2) = G{xi, Ai, Bi; X2, A2, B2) ± G{xi, Ai, Bi, X2, -A2, -B2). 

(22) 

In the subsequent sections we will analyze the properties of G. 

To close this section let us point out that the propagator G can be associated 
with a field propagator 

G(a:i,pi;a:2,P2) = (0|r{V^i„B, (^i)V'a.U(^2)}|0). (23) 

Similarly, Green functions for any number (p) of operators ■0a,b(^)) ^^<i oP" 
erators (pa{y), can be defined (with the restriction rii — ±1) according to: 

{o\T{rA\,BA^i) ■ • •^A;B,(^p)0a.(?/i) . ■■<t>aM)m = 

lim/p[t/]0„,(yi)...0„,(y,)e-^t'Jl+P'^% (24) 
together with the boundary conditions, 

AiXiBj ff rij = -1 

' (25) 

Aii-XiYBj a m = +1 



^ We define A > 0, A e SU(2) when A = qq + ia.f, (eq + = 1) is written with 
ao > 0. 
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4 4 



X, = -—^ + i^^Ta e SU(2), 
for small spheres Si around Xi (Cg = — (37r^/e^) log(e)). 

If n > 1, the operator ij/\x) could be defined as a suitable limit of ip'^{yi) . . . ip^{yn) 
when Ui x. However we will not pursue this point here but instead, from 
now on we will consider only ■0='=^ fields which will be denoted as ■0+^ = ■0, 



4 Symmetries of the propagator 



The symmetries of the propagator can be found by performing different changes 
of variables, within the path integral, which are invariances of the action. In 
this way we learn how this symmetries act on the internal variables Ai, Bi. 

Translations The change of variable ?7(x) U{x + a) does not affect the 
action but affects the boundary conditions. Trivially, it follows that 

G{xi + a, ^1, Si; X2 + a, A^, B2) = ^(xi, ^1, B^; x^, ^2, B^). (26) 

Isorotations With the change of variables IJ = AUA'^ , {A e SU(2)) we 
obtain 



S[AUA^ 
AlA^ 



U 



Si 



S[U], 
1 at 00, 

A U\s3A'^ = AAiXiBIA'^ = {AAi)Xi{ABi)\ 

A U\ssA^ = ^^(-X2)^S|At = (AA2)(-X2)t(AS2)t(27) 



Thus 



G{xi, AAi, AB,; x^, AA^, AB^) = G{x^, A, B^; ^2, A2, B^). (28) 

Lorentz transformations To describe particles of half-integer spin we 
must consider the universal covering of the Euclidean Lorentz group, 
namely SU(2) x SU(2). To each pair of matrices 81,82 G SU(2) corre- 
sponds a matrix A G SO (4). The transformation x ^ x' — Ax can be 
written in matrix form 



X — > X' — 8\X82i X 



+ i^TaeSU(2). 



\x\ 



Now let us perform the change of variables 

U{x) U{x) = U{A-^x), 



(29) 



(30) 
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which leaves the action invariant. The boundary conditions transform as 

f/(X')|^,3 = UiX)\,s = A^X^Bl = A^SlX[S2Bl = iA^Sl)X^iB^Sly , 
U{X')l,^ = U{X)\,, = A,{-X^,)Bl = {A,Sl){~x',^){B,Sl)\ (31) 

which gives 

G{Axu A^Sl B,Sl; Ax^, A,Sl B^Sl) = G{x^, A^,B,; X2, A^, B^). (32) 

Parity The pion is pseudo-scalar, so under parity U transforms as 

U{x\ x) U{x\ x) = U^{x\ -x). (33) 

As a spatial reflection maps Xi 2 XI2, a parity transformation simply 
interchanges A and B in the boundary conditions. So, parity invariance 
of the action implies that 

G{xl, Xi,Ai,Bi; X2, X2, A2, B2) = G{xl, -xi, Bi,Ai; x\, -X2, B2, A2). 

(34) 

Let us proceed now to find how these symmetries constraint the form of the 
propagator. Because of translational invariance G depends on the difference 
Xi — X2- The vector xi — X2 can be split into a radial and an angular part as 
follows 

{\xi — X2I 

G is a smooth function of X12, Ai, Bi, A2, B2. Any such function can be de- 
composed as a linear combination of the matrix elements D-'^^,. This can be 
seen from group theory or because the matrices are the spherical harmonics 
on S^. Thus, we expand G as 

G{\xi-X2\,Xi2,Ai,Bi,A2,B2)= Yl vUi:(^i,(^i',\x2 - 

ji, Wi\,Wi\<ji 

K.', (^2)/^^:., {B2)D^:,^, {X,2). (36) 

where j runs over \Z>q, i.e. integer and half-integer representations. From 
definition (^) we see that 

G(xi, Ai, Bi- X2, A2, B2) = G{xi, -Ai, -Bi; X2, -A2, -B2) 
where the implication follows from the formula D^[—A) = {—)'^^D^{A). 
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Performing independent variations of A, Si, S2 in eqs. (|28|), (|32D and using 
standard arguments about the coupling of SU(2) representations one finds 
that the rj indices (cri(T2cr3(j4) are coupled to zero (isorotations), as well as 
the indices (o"2, cts, cts) and (cr^, CT4, — CT5) (Lorentz transformations). Thus, the 
most general G satisfying all these symmetries is 



-X2|,Xi2, A, 51,^2,52) = 

(2-^ + l)^<^<ai<^a-40-^C'm2mi^jij2i3j4i5/(k2 — 

all <j,<j' ,j,m,I 

I 




J3 32 J5 



V^3 



' a' a' 



3i 32 I I [ is 34 I 
(Ti 0-2 mi / I 0-3 0-4 m2 



(38) 



where use was made of the Wigner's 3j symbols (which couple to zero three an- 
gular momenta) and of the matrices Cmm' = i.—V^^'^m.-m'- Useful properties 
of C are C = C*, C-^ = C* = {-f^C. 

Parity conservation imphes Vjij2j3j4j5i = Vj2jij4,j3j5i (cf- eq. (^If)). 

The transformation laws of this section can also be defined over the fields 
4'A,B,'4'A,B by means of unitary operators U satisfying: 

Isorotations: 

Ua^Ai,Bi(x)U^ = IpAAuABii^), 

Ua^Ai,Bi(x)U^ = iJAAuABAx). (39) 

Lorentz transformations: 

UAV'Ai,Bi(x)Ui = ^^^gt^B^gt(Ax), 

UA^A2,B2(x)Ui = ^A^gt B^gt(Ax), (40) 



where A = {Si, 6*2) . 



These transformations suggest changing to a basis of fields which transform 
under a finite dimensional representation of the Lorentz and isospin groups. 
Finite dimensional irreducible representations of the SU(2) x SU(2) group are 
labeled by two numbers ji,j2 ^ |^>o- If A = (Si, 5*2) the representation 
(j,0) is given by V^{A) = D^lSi) and (0,j) by &{A) = D^{S2). The matri- 
ces are the usual rotation matrices representing SU(2) and so, unlike the 
Minkowski case, this representations are unitary but not-faithful. Finally the 
{jiih) representations are given by (ji,0) ® (0,j2)- The fields transforming 
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under these representations are given by: 



/ 



Ji h I 
(j'l ^ 



(41) 



A.B 



and the same formula for if). The integration over SU(2) is performed using 
the group invariant measure and Zj^j^j is a normahzation factor to be chosen 



later. The transformations ([39|), pO|) now read: 



/Jlj2l 



'CTiCT2mi 



X 



m'jTnj 



\ 'Tcr-icr2inj 



X), 

,(x), 



tTi(T2mi > 

jlj2l 
(Ti(T2 mi 



X 



)U 



(42) 



which means that V'^ii^am/ (^) transforms under the Lorentz group representa- 
tion (ji, ^2) and 4'i\^a^mii^) under the representation (72, ji)- In addition, both 
fields have isospin /. 



and definition (O) we can obtain the general form of the prop- 



Using eq. 
agator 



E 





32 




'1. 


J4 


J6 1 




0-2 






< 





(43) 



where 77 can be written in terms of 7] and the 6j symbols []TH| as 



Vjlj2j3jd6l 




(44) 



The spins associated with these fields are obtained by composing ji and 72 
and so, in principle, they range from |j2 — ji| to ji -|- j2- The propagated spins 
are determined by the propagator in ([l6|) . In the next section we will show 
that (|TB|) has a new symmetry which leads to the propagation of particles with 
spin equal to the isospin /. 
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5 Skyrmion propagator 



Now, we will approximately evaluate the propagator in the large \x2 — Xi\ 
limit. Performing a Lorentz transformation, the 4- axis can always be taken 
along X2 — Xi, i.e. X12 = 1. By means of eq. (^31) we can extend this result to 
arbitrary values of X12. 

Having fixed X12, we will analize what is expected for different values of 2 
and Bi 2- 

li Ai = Bi = A2 = B2 = I, the situation is that of section |^ where we 
evaluated (numerically) the saddle point. This saddle point will be the base 
for our analysis and will be denoted as Uo{t,x). One should remember that it 
behaves as —Xl around X2 and as Xi around xi. Note also that Uo{t,x) (as a 
function of x), in the time interval between X2 and xf, is almost equal to the 
static soliton profile, which will be denoted as Usk{x). 

Suppose now that A2 = B2 = A, Ai = Bi = 1. This corresponds to performing 
an isorotation in the vicinity of X2- We know that if ?7 is a solution to the 
equation of motion, then, AU A^ (with constant A) is also a solution, that is, 
the function AUqA^ satisfies the equations of motion around X2- It also satisfies 
the boundary condition on 9-B|. Moreover, this behavior can be matched with 
the behavior around Xi by means of a time dependent ansatz of the form 

U{t,x) = A{t)UoA\t), A{xt) = A, A{xl) = 1, (45) 

which satisfies f/ = 1 at infinity; upon substitution in the action we obtain 
an equation for A{t). The corresponding solution will represent, through eq. 
(H3), the creation at X2 of a skyrmion isorotated in A, with profile AUsk^^ , 



which slowly rotates until it is annihilated unrotated at Xi. This solution is 
expected to give a good approximation to the saddle point. More generally, for 
A2 = B2, Ai = Bi the saddle point will be approximately given by minimizing 
the action with respect to the ansatz (45) with ^4(0:2) = A2, A{xf) = Ai. This 



is by now a standard calculation in skyrmion physics which gives free motion 
over SU(2) with a moment of inertia 9", obtained from Usk{x) A better 
approximation to the field propagator can be done by integrating over all 
possible A(t) giving the free propagator over SU(2), from A2 to Ai, which 
corresponds to quantizing the isorotational (or rotational) zero-modes. 

Now, we consider the case A2 B2, Ai = Bi = 1. In this case we also have that 
A2U0B2 is a solution in the vicinity of X2 that satisfies the correct boundary 
condition at (9-B| but here, a time dependent ansatz does not suffices; for such 
an ansatz, the behavior A2U0B2 around the vicinity of X2 would be extended 
to infinity on the hyperp lanes x^ = const, x^to a. value A2B2 7^ 1 that does 
not match the required asymptotic behavior of the fields. Recall that, around 
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X2, Uo is unlocalized for |a; — a;2| < fi, while it is localized on the North pole for 
|x — X2I > Ti (see the discussion at the end of section ^). Then, if we start with 
a function such that, very close to — a:2| = e, behaves as the solution A2UQB2 
and this behavior persists up to the localization radius, the topological density 
will not be localized on the North pole of the sphere |x — X2I = r;, but on the 
direction x — X2 satisfying X2 = B2A2. This behavior represents a soliton that 
starts moving on a direction other than towards Xi. Therefore, the matching 
with the behavior around Xi where the soliton is required to be annihilated will 
lead to a function with action larger than the action associated with a function 
representing a soliton that starts moving towards Xi, that is, we expect that 
the configuration which minimizes the action behaves, in the unlocalization 
region around X2, as 

U{x) = A{t)UoB{t)\ A{e) = A2, B{e) = B2, A{ri) = B{ri) = C2 (46) 

where r = |x — X2I, e is the radius of the small sphere we extracted around 
X2 and ri ~ (c/tt)"^. Note that this behavior, at r = r;, correctly localizes the 
topological density on the North pole, though isorotated by C2. The value of 
C2 must be chosen so as to minimize the generated action in the unlocalization 
region. Loosely speaking, (C2, C2) will be the closest rotation to the Lorentz 
transformation (742,-82)- A more quantitative answer is obtained by inserting 
the above ansatz in the action, leading to 



S = A-n'^ J dr h(r)Tr (AA^ + BB^ 



e<T<ri 

where the dots represent derivation with respect to r, and we made the approx- 
imation Uq ~ —Xl- After a change of variables u = u{t) such that ii = l/hi^r), 
which cancels the factor /i(r), the expression (|47|) represents free motion over 
SU(2) X SU(2) for {A{u),B{u)). For small r, the solution to -it = l//i(r) is 

At M ~ 'u(r/) we must have A{u) = B{u) = C2, then the action is equal to 

u{ri) - u{e) 

SAB = 2Arccos(iTr(AB^)). (48) 

Where sab is the distance between A and B measured with the standard 
metric of SU(2). For A2, B2 fixed, a calculation reveals that this action is 



minimized by C2 = B2y B2A2. Note that in SU(2), the square root is well 
defined up to a signQ. This part of the saddle point action contributes to the 

^ except for the square root of minus the identity since any matrix A = ia.f, a.a = 1, 
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propagator with a factor exp{—{s\^Q^ + s\^q^) / {u{ri) — ^(e)). In a general 
situation, where Ai ^ Bi, the same analysis is also valid around Xi. The 



results are best expressed in terms of the minimizing variables Ci = Biy bJA^ 
and the variables Di, defined by the relations 

Ai = CiD,, Bi = aDl (49) 



which yields Di = y BjAi. For the minimizing Cj's we have 

4a + sic. = = 2Arccos(iTr(Df)). (50) 



Now, the behaviors (at r^) around Xi and X2 can be matched by means of a 
time dependent ansatz whose minimization will complete the construction of 
the saddle point and will represent a slowly rotating soliton from C2 to Ci. 
Again, a better approximation is obtained by considering free propagation 
over SU(2) from C2 to Ci (whose expression can be found in [0); including 
the factors generated in the unlocalization region, we obtain the following 
expression for G 



G{x,,Ar,B,;x2,A2,B2) = (7i(TrDi)g2(TrD2) ^ Dkm(Ci)DLk(C5) 

I 

-Mci|x2~xi| 

1x2 - Xi 



e-*l---l^ (51) 

V M I 

where we also include the factor exp —Mc\\x2 — Xi\ coming from the classical 
action, and the factor proportional to \x2 — resulting from the inte- 

gration over the translational zero-modes. The function gi only depends on 
Tr(Di), since for an SU(2) matrix Tr(D^) = (Tr(D))^ — 2, and contains the 
gaussian factor exp{-{s%c, + ■s|,cJ/(«(^/) " ^(e)) (cf- eq. (|0D). 

Finally, to obtain the propagator G, we note that A —A, B —B is 
equivalent to C -> -C, D ^ D, and D^{-C) = {-f^D^{Cy, therefore, if we 
define G in eq. ([2^ ) using the plus (minus) sign the only surviving terms are 
those with integer (half-integer) /. We will show below that / is the spin and 
isospin, thus choosing the plus (minus) sign will lead to integer (half-integer) 
spin (equal to isospin) particles. 

Note that G in equation (^TJ) enjoys the two following properties: (1) G satisfies 
the parity requirement ( P^ ) (invariance under Ai Bi) since this is equivalent 
to Di ^ D\ and for an SU(2) matrix Tr(D) = Tr(D'''), (2) G is invariant under 
the new symmetry Di — > RiDiR-. 

satisfies A"^ = —1. 
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Now we will show that the propagator of eq. ([5T|) propagates spin equal to 
isospin particles. 

To find the propagator of the fields ip^^^'-^ we use eq. (^If) to obtain: 





k 
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ni'j j 




4 


nil j 



Crr' m C rrL rro C rrL rr-j C rr' 



Al 2^1,2 



^i(TrDi)g2(TrD2)DL(Ci)DLk(C5 



g— Mj|x2-xi| 



[ 2ZI 
V M 



X2 - Xi 



3/2 ■ 



(52) 



where we introduced Mj = Md + /(/ + l)/(253). Let us first perform the 
integral over Ai, Bi. It is convenient to write everything in terms of Ci and 
Di, the factors involved are: 



Ji J2 / 



Co-; CTi e'er' a-2 ^ 



X 



,(CiZ^O^iL'(Ci^I)^?i(TrDODL(Ci). 



(53) 



To change variables we must evaluate a jacobian which can be calculated 
parametrizing the matrices as, A = + ia.r with (oq + a? = 1) and similar 
expressions for B, C and D. The invariant measure is written as 



dA 



dai da2 das 
ao 



(54) 



The jacobian is the determinant of a 6 x 6 matrix which upon evaluation gives: 

J dAdB = J dCdD2(TrD)l (55) 

We absorb this factor in gi defining gi = 2(71 (TrDi)^. Since cji is a class 
function of Di, that is gi{SDiS'^) = gi{Di), it can be expanded according to 



(56) 



We are now ready to evaluate the integrals since the integral of the product 
of three D matrices is given by the 3j symbols ||15[. The result is 
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Jl J2 I \ Jl J2 I 



o'l 4 mi j \ a'l a'i k 




( , 



Jl 32 k 1 / jl J2 ^1 



.(57) 



The two first 3j symbols give dijdmjk/ {'^I + 1), the three remaining symbols 
can be simplified by means of the 6j symbol to give 



gih ^ ^ r \ 3i 32 I \ { Jl J2 I 



(58) 



The integral over ^2, B2 can be performed in a similar way. Putting all the 
pieces together we obtain: 



1 Cm'^rni^II' G 



-Mi|x2-xi 



jl J2 I \ h ji I 



Gi 02 ra \ (73 (74 m 




{jm\x2 - Xi 

jl ^2^11 ^3 j4 / 

jl 32 k js j4 ^2 



3/2 



(59) 



Choosing in this equation 

y 

and the same for Zj^j^j/ (replacing gn-^ by (72/2) arrive at 
/ 




(60) 



Jl J2 / W J3 31 I \ e-^'l"^"" ^l 

^m'm^in'jTnj^II' 



ycTi a2 m J y 0-3 0-4 



3/2 ■ 



(61) 



Finally, from eq. (^3[) we can read the value of 77 and then extend the propa- 
gator from X12 = 1 to an arbitrary X12, obtaining 



(o|T{^^^/:,e/^/„j|o) = C,„'^C„^™,5,,.I?^;,,(Xi2)I^^:.,(Xi2 

/ 




j3 34 I 



Mi|x2-xi| 



CX3 m'y (^^|a;2 -xil) 



3/2 



(62) 



This expression coincides with the free field propagator we evaluate in ap- 
pendix ^ namely, formula ( |A.15| ), where the index / in (p^ ) replaces the 
index j in that formula. The index j in (|A.15|) is the spin of the one particle 
states we use in the appendix to construct the free fields, on the other hand. 
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I is the isospin of our fields (cf. eq. implying that the propagator G de- 

scribes a tower of spin equal to isospin particles (with mass Mi). Furthermore, 
the Green functions obtained in the appendix correspond to free fields that for 
integer (half-integer) j commute (anticommute) at space-like separations, that 
is, the fields we define by means of the Green functions ( |23D and eq. (^) obey 
the spin-statistics theorem. For example, one can calculate from eq. ( |62D the 
propagator in the (1/2,0) © (0,1/2) representation which coincides, at large 
distances, with the usual Dirac propagator for a spin 1/2 particle. 



6 S-matrix definition 



In order to compute scattering amplitudes one should be able to extract the 
S-matrix from the Green functions previously defined. As a first step we apply 
definition (^Tj) to obtain Green functions for the fields ip^'^^^ and t/'-'i-'^ (remem- 
ber that if) and ip refer to singularities of topological number -|-1 and —1 re- 
spectively). Then, the Green functions are analytically continued to Minkowski 
space-time by the replacement — — ix". The large distance propagator, ob- 
tained in the previous section, allows us to identify the asymptotic states to be 
scattered, in the non-trivial topological sector, as composed by spin equal to 
isospin particles. Now, all we need is the corresponding LSZ reduction formula, 
which is obtained in the appendix and reads 

out(pia-iSi, ■ ■ ■ ,p„a„s„|gi(TiSi, ■ ■ ■,qicrisi)in = disc, parts + (iZ"^/^)''+^ 
J dVi ■ ■ ■ d\u^^^'^,^ (pi, ai) • ■ ■ u;|^;^, (qi, exp |^i J2 PkYk - i ^^Xr^ 
+ml)--- (□., + mf)(0|rV^^;L';,^(yO . ■ -CS JxOlO), (63) 

where cxj stands for the spin projection z in the rest frame and s is the isospin 
grd projection. The wave functions u, u are defined in the appendix; when 
considering the scattering of antiparticles, u, u must be replaced by f , v. The 
Green function in eq. (^) is obtained from ( p^ ) and (|4T|). An intuitive under- 
standing of this Green function can be achieved using the line defects picture 
of ref. 1^. There, the leading contribution to the path integral ( P^ is identified 
with configurations described by open line defects emerging from xi, x„ plus 
a gas of closed defect tubes. If one brings, e.g., xi far from the other points, 
the sum over the open line defects associated with xi will give a skyrmion 
free propagator, from Xi towards the interaction region. This propagator has 
the pole (in p^) required for a non-vanishing scattering amplitude and also 
contains all information about the boundary condition around Xi. 

As it is noticed in the appendix, the different fields that can be built from 
the same particle states are not independent and one should choose only one 
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iphjii represent each kind of particle. For example the A (j = 3/2) can 
be represented by a (3/2,0) field or by a (1, 1/2), as in the Rarita-Schwinger 
representation. Choosing one or the other amounts to a field redefinition and 
hence the S-matrix should be independent of this choice. This assertion can be 
checked, within the saddle point approximation, taking into account that all 
information about the boundary conditions is contained in the free propagators 
( p^ ) attached to the points described in the previous paragraph). The 

explicit dependence on (ji, ^2) is cancelled by the choice of normalization factor 

It is interesting to note that, for a (1/2,0) field, the reduction formula we 
give contains an operator (□ +m^), instead of the usual (i^ + m). This comes 
about because we are using a two-component spinor, see Jl^ for a discussion 
on this point. For a two-component spinor we need ip^^'^ and dti^^^"^ ° to 
recover the operators a and 6^ In the case of spin 1/2 fermions, the usual 
alternative approach is considering two fields, namely ip^^'^ ° and tp^ (or in 
general ip^^^'^ and ip^^^^). Using these two fields (4-component spinors) one can 
obtain a and with no need of the field time derivative. Then, in that case, 
the LSZ reduction formula displays a first order operator. The same can be 
done for higher spin fields but the operator appearing in the LSZ formula will 
be of higher order and more cumbersome to use than the operator (□ -|- m^) 
associated with a single field representation. 



The Green functions we defined in eq. (|2^) imply that our pion field is not 
just the fiuctuating part but the sum of the classical configuration plus the 
fiuctuating part. In fact, this is the natural thing to do, the separation be- 
tween classical and fiuctuating parts is just a consequence of the saddle point 
approximation, which may not always be adequate. 

Precisely, the suggestion that the pion field be identified with the classical 
soliton plus fiuctuations, instead of identifying it just with the fiuctuations, 
led to the solution to the Yukawa problem (see refs. |TB[ and ]I^). 



The Yukawa coupling is obtained from the vertex 

(O|T{<',^i,(xi)e;,^;^li(^2)0a(yi)}|O^ (64) 

At large distances (|xi — X2I large) the relevant configurations are described in 
terms of open line defects which can be interpreted as skyrmion world lines. 
Then, in this case, the saddle point calculation is similar to that of reference 
where the correct Yukawa coupling was obtained resorting to a rotationally 
improved ansatz. In our formalism, the calculation would be anologous to the 
evaluation of the propagator performed in section |^. 

At short distances (|xi — X2I small) the saddle point no longer resembles a 
static solution in the interval between creation and destruction of the soliton, 
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so one should expect differences with the approach of . However we will not 
pursue this point here since at short distances the Skyrme lagrangian is not 
adequate. However, this discussion may be relevant in connection with other 
solitonic models. 



7 Conclusions 

In ref. Frohlich and Marchetti proposed a quantum field theory for skyrmi- 
ons, defining Euclidean Green functions in terms of path integrals over singular 
fields. 

In the present paper we show that by fixing the boundary conditions around 
the singularities it is possible to study the physical content of the Skyrme 
model and their dynamical implications in a manifestly covariant way. In par- 
ticular, we obtain the covariant baryon propagator. This treatment may also 
be relevant when studying other solitonic models as for example monopoles. 

Our method is based on defining Green functions subject to a minimal set of 
boundary conditions that close under four dimensional rotations and isospin 
transformations. The fields ipA,B defined in this way are labeled by two SU(2) 
matrices and carry an infinite (reducible) representation of the isospin group 
and four dimensional rotations. Equivalently, we can define a set of fields 
^jd2i transforming under finite dimensional irreducible representations. Here, 
the spin of the field ip^^^^^ ranges, in principle, from \j2 — ji\ to ji + j2 (/ is 
the isospin); for these fields we define the corresponding Green functions and 
derive the LSZ formulae to compute the S-matrix, thus mapping a complete 
relativistic field theory. The S-matrix constraints of locality and unitarity are 
expected to be satisfied. An heuristic argument could be given along the lines 
of ref. using the Osterwalder-Schrader reconstruction theorem. 

The numerical calculation of the saddle point shows that near the singularities 
we have a region where the chiral symmetry is unbroken, while it is broken 
away from the singularities. These regions correspond to the unlocalization 
region, where the topological density is isotropically distributed in Euclidean 
space-time, and the localization region, where the topological density gets 
concentrated, signaling the presence of a soliton. 

In order to compute the path integrals, we show that, in the localization region, 
and when the skyrmion is created at rest, the relevant paths can be character- 
ized by (non- covariant) variables C = B\J obtained by minimizing the 
action generated in the unlocalization region. 

The covariant propagator is evaluated by using symmetry arguments and inte- 
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grating over zero modes in the localization region. The result obtained shows 
a new symmetry under the transformation D RDR\ where D = V B'^A. 
This symmetry is seen to be responsible, upon an explicit calculation, for the 
spin equal to isospin skyrmion spectrum of the model. 

Another consequence that shows up, upon fixing the boundary conditions, is 
the introduction of additional topology, which leads to consider the extension 
of a map U : x ^ SU(2) to a map U : x D2 ^ SU(2) to study the 
different possible inequivalent quantizations of the SU(2)-Skyrme model 0. 
There are two possibilities, one of them corresponding to the baryonic spec- 
trum composed by half-integer spin (equal to isospin) particles, which are seen 
to satisfy the spin-statistics theorem. 

Finally, we think that our method also provides a formal base for the work in 
ref. where a skyrmion field theory is advocated (together with an improved 
approximation scheme) to solve the Yukawa problem, that is, reproducing the 
pseudovector pion-baryon coupling from skyrmion physics. 
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A 



In this appendix we summarize a method due to Weinberg to define 
arbitrary spin field propagators in Minkowski space. 

We start by considering states containing one particle (with mass m) at rest, 
having spin j and spin projection a: |(m, 0); jcr) and then we define the states 

b;ja>= ./-[/(Lp)|(m,0);ja), (A.l) 
V Po 

where U (Lp) is a boost that maps (m, 0) into p, which is defined by the SL(2, C) 
matrix 

Lp = expiep.r). (A.2) 

where p is the unit vector p/\p\, and sinh(^^) = \p\/m. Note that LpLp = 
(j>o + p.T)/m. The Lorentz transformations act over the states ( [A.!]) according 
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to: 



UK\p;3a) = J-U{AL,)\{m,Oy,ja) = J -U{L^p)U{L-^lAL,)\{m,Oy, ja) 
Po V ^'o 



^^^^'Di,ALAlAL,)\Ap;ja'). (A.3) 



Po 



The key point is that L^^ALp is a rotation and rotations act on the states 
|(m, 0); jcr) by means of the matrices, as usual. 

Particles in the state |p; jcr) are created from the vacuum by operators aj^^ 
which, in view of eq. ( [A.3|) , satisfy: 



V Po 

UWp.M = \^^{CD^{L-'A-'L^p)C-'U,al^^,. (A.4) 

V Po 

Using the operators a, a^, and the corresponding operators b,b^ for the an- 
tiparticles, fields in any finite dimensional representation (ji, ^2) of the Lorentz 
group (with |j2 — ji\ < j < \ji +J2I) can be constructed through the formulae 



Cera' 



Jl J2 J 



The sign (— )^-'^ is needed for the ip and ip"^ to commute (or anticommute) 
at space-like separations. At time-like separations they do not commute since 
they are built with the same operators a, a^. This means that the various 
relativistic fields we can build from the same spin j particle states are not 
independent; we must choose one of them to represent these particles. The 
field transforms as 

UACX(x)Ui = 2^i^,^;(A-^)P^^^^, (A-i)<f^, (Ax) (A.6) 

while ip obeys a similar relation but with ji,j2 interchanged, as corresponds 
to the adjoint operator. Now we can readily evaluate the propagator 



(0|T{^(x)^^f„^V^(2/)^^f„J|0) 



H'HOmyyr^^rnM^y^Tn^JO) VO > ^0 

(A.7) 
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giving 



{2n 



^ 3i 32 j 
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J3 J4 J 



^^(^0 - yo)e-'P(^-^) + (-)('j^+'j^)^^(yo - xo)e'P(^-^)l (A.8) 



For a boost &{Lp) = V^{Lp^); then, using the property of the 3j symbols: 



Jl 32 33 



\ m^^ m2 I 
This expression can be transformed into 

(27r)3^cj. "^"^i™'! (-^P-^p) -^^4™^ {LpLpjCa'a 



Jl J2 J3 
777-1 "^2 ^3 



(A.9) 




J3 J4 J 



7713 "^4 



e(a;o - yo)e-^P(^-^) + {~T'''^^'''^9{y^ - xo)e'P(^-y)l . (A.IO) 



Taking into account that 



V^I\LpLp 



LpLp = —(po+P-t) 

777 



(A.ll) 



and that the elements of the matrix are linear combinations of the elements 
of P^/^®. . .'^V^^'^, {2j times), it follows that is a polynomial inp^ of order 
2j. In the integral (|AlOD we can "take out" the polynomials V^{{pQ+p.f)/m), 
if we replace by — i9^. So the final expression for the Minkowski propagator 
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ido + T.d))Apix-y), 



(A.12) 



where we introduced the scalar Feynman propagator A^. When extracting the 
time derivatives we did not applied them to the 6 functions that define the time 
ordering. Thus, the above propagator is the so called covariant propagator. 
The non covariant pieces should be automatically cancelled by non covariant 
vertices in the interaction Hamiltonian when using the canonical formalism 
(see [0). In our case, since the path integral used to define the propagator (eq. 
([T6|)) is Lorentz invariant, we also expect to obtain this covariant propagator. 
As the {—ia^dmU.) are polynomials in 5^, this expression is suitable for 
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analytic continuation into Euclidean space. Defining the Euclidean propagator 
as Ae{x4,x) = Ap{ix4,x) we obtain the analytic continuation 



+ i^^-^^))^m4mi(-;^;(^o + i^^-^"'))^E(x - y). (A.13) 
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For large distances Ae{x) behaves as 
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and to leading order in l/(m|x|) we can replace 9^ by — mx^/x, giving 
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where we introduced the matrix X = X4 + ix.r. 



Another useful calculation that can be done with the free fields (|A.5|) is the LSZ 
reduction formula. The procedure is well-known [^. First we must express a, 
a\ b, in terms of the fields ip. This is an easy task which gives, for example, 
the operator as: 



,t(m) 
''pa 




(2^)3/2^ ^ MS-^ - i^t^iiiLJ • (A.16) 



The LSZ reduction proceeds now by using the asymptotic condition for the 
interacting field 

(/|^(t)|z)*^°^Z-V2(/|^(in)|,), (A.17) 

to replace if)'^™'^ by ip in the formula for a^^. The renormalization factor Z will 
absorb the factor exp(— C^) introduced in eq. (pIBI) . Then, using the relation 
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t — >— 00 J t — ^+00 

—00 

we arrive at the reduction formula. The result for as well as for the other 
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operators is 
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(out|b(7*)|in) = (out|bJ°)|in) +iZ 
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= disc. + iZV2| d^xe'P^vjif^^(pa)(n + m2)(out|7/;j^f^^ 

where for brevity we defined -u, -u, f, -u which can be interpreted as wave func- 
tions of incoming and outgoing particles. These formulae allow us to extract 
the scattering amplitudes from the Green functions we defined. 
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